Random-Mass Dirac Fermions in an Imaginary Vector Potential (II): 
Long-Range Correlated Random Mass 
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In the previous paper, we studied the random- mass Dirac fermion in one dimension by using the 
transfer-matrix methods and by introducing an imaginary vector potential in order to calculate the 
localization lengths. Especially we considered effects of the nonlocal but short-range correlations 
of the random mass. In this paper, we shall study effects of the long-range correlations of the 
random mass especially on the derealization transition. The results depend on how randomness is 
introduced in the Dirac mass. 
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PACS: 72.15.Rn, 73.20.Jc, 72.10.Bg 
I. INTRODUCTION 

In the previous papers we studied the random- 

mass Dirac fermions in one dimension by using the 
transfer-matrix methods (TMM) and imaginary-vector 
potential method (IVPM). We calculated the density of 
states, the typical and mean localization lengths, and the 
multifractal scalings as a function of the energy and cor- 
relation length of the random mass. The results are in 
good agreement with the available analytical calculations 

. We also obtained the relation between the correlation 
length of the random mass and the typical localization 
length for the short-range correlations 

In this paper, we shall consider effects of the long-range 
correlations, especially on the delocalization transition. 
This problem is interesting. Very recently finite mobil- 
ity edges are observed in the one-dimensional Anderson 
model with the long-range correlated potentials though 
it is widely believed that (almost) all states are local- 
ized in spatial one dimension Similar phenomena 
are observed also in the aperiodic Kronig-Penney model 
In contrast to the above models, the system of the 
random-mass Dirac fermions belongs to the universality 
class of the chiral orthogonal ensemble and therefore the 
extended states exist in the band center though the oth- 
ers are localized for the white-noise random mass. Long- 
range correlations of the disorder may change this fea- 
ture. This problem has not been addressed in detail so 
far. As the model is closely related with the random-bond 
XY model, the random Ising model, etc., the results are 
interesting both theoretically and experimentally. 

This paper is organized as follows. The model and nu- 
merical methods are explained in the previous paper . 
In Sec. 2, we shall explain how to make a long-range corre- 
lated random Dirac mass numerically. We shall consider 
two types of the telegraphic random mass, i.e., in the 
first one magnitudes of the random mass are long-range 
correlated random variables (LRCRV) with fixed inter- 



val distances between kinks whereas in the second one 
interval distances between kinks are LRCRV with a fixed 
random mass magnitude. As we show in later, physical 
results are different in the above two cases. As far as we 
know, this possibility has not been found in the previ- 
ous studies on the random systems. In Sec. 3, numerical 
results are given. We calculate the localization lengths 
directly by using the TMM and IVPM though in most 
of studies of one-dimensional random models elaborated 
techniques like the Hamiltonian mapping, the renormal- 
ization group, etc. are used for calculating the Lyapunov 
exponent. Before obtaining the localization lengths in 
the present methods, we have various landscapes of the 
wave functions themselves and obtain a physical picture 
of the "metal-insulator" transition from them. Section 5 
is devoted to conclusion. 



II. MODELS AND LONG-RANGE CORRELATED 
RANDOM MASS 



Hamiltonian of the random- mass Dirac fermion is given 



by 



n 



h = —ia^dx + m{x)a'^ , 



(1) 
(2) 



where a are the Pauli matrices and m{x) is the tele- 
graphic random mass. In this paper we shall consider the 
long-range correlated configurations of m{x) and study 
delocahzed states by the TMM and IVPM. Numerical 
methods which generate long-range correlated potentials 
was recently invented by Izrailev et. al. and Herbut [^,0 . 
We briefiy review it for we shall use it for the present 
studies. 

Let us consider a spatial lattice and random potential 
e{n) sitting on the sites which has correlation x("-)7 



[ e(m) e(n) ]cns = C x(|m - n\), 



(3) 



1 



where n and m are site indices, and C is a constant. Here 
we suppose that xC'^) is given and define c(/c) from x(n) 
as follows, 

m^m'^', (4) 

where x{k) is the discrete Fourier transform of x(n); 



(5) 



We also introduce another random site potential r{n) 
which is distributed uniformly in the range [—1, 1] with 
the delta-function whitc-noisc correlation, 



r{m) r{n) ]o„s oc (5,„,„. 



(6) 



Then we can construct random potential e(n) as follows 
by using the white- noise potential r(m) and the inverse 
Fourier transform of c{k), 



e{n) = C r{m) c{m — n), 



(7) 



way that the typical value of e(m) does not change. Then 
as a result, the strength of correlation does not change 
for various system size L if we fix |m — n\ in Eq.(lO). 

The above method generates the random potential at 
each site (Fig.l). Therefore, if we fix the distances be- 
tween kinks and vary the value of m{x) for generating 
random m(x), we can directly use this method. We 
can extend TMM easily to the case of random values of 
|m(x)|. Numerical studies of this system will be reported 
in Sec. 3. A. 




where C is another constant. However, we cannot per- 
form the infinite summation in the numerical calculation, 
and therefore we replace it with the finite summation as 
follows. 



k/2 

e{n) = C r(m) c{m — n). 

m=-k/2 



(8) 



where A: is a finite number which is much larger than the 
number of sites.0 We set k ~ 10^ -I- (number of sites) in 
our numerical study. 

We can easily verify that the correlator (^) is satisfied 
from these variables. If we choose c(m) as 



c(m) 



TO 



(9) 



then the correlator of the random potential becomes 

1 



e(m) e(n) 



\2a-l ■ 



(10) 



With this procedure, we can generate power-law corre- 
lated random potential. 

In the actual calculation, there is a subtle point in the 
normalization of the random potential e(n) in Eq.(p^) 
when we vary the system size L We determine the 
normalization of e{in) (namely the value of C ) in such a 



*We use the periodic boundary condition for the numerical 
calculation in subsequent sections. 



FIG. 1. An example of m{x) 



On the other hand, we can also use the above method 
for generating random distances between kinks with fixed 
value of \m{x)\. This is the subject in Sec.S.B. Physical 
results are different in the above two cases. 



III. LONG-RANGE CORRELATED DISORDERS 
A. Random mass with random magnitude 

We focus on the system with power-law correlated ran- 
dom mass in this subsection. As we explained at the end 
of Sec. 2, we first consider the random mass of random 
magnitude with fixed distances between kinks. Prototype 
of the configuration is given in Fig. 2. In this case m{x) 
is generated rather similarly to e(n) in Sec. 2 though the 
continuum space is considered here instead of the spatial 
lattice. As in the previous paper we use the IVPM 
and TMM in order to calculate the localization lengths. 

We first calculate the correlation [m(x)m(0)]ens numer- 
ically which is expected to exhibit the power-law decay. 
The result is shown in Fig. 3. The correlator actually 
shows the algebraic decay, though the decay powers are 
slightly smaller than the expected values 2a — 1 (denoted 
by Oipw)- We think that this is due to the finiteness of 
numbers of kinks and/or finite Fourier transformation in 
the scheme which we explained in the previous section. 
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FIG. 2. An example of m{x) in the case of power-law cor- 



related random mass : We set a„ 



: 0.5 and 2000 kinks. 



mobility edge at a finite energy which roughly equals to 
Max {\m{x)\}. 
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FIG. 3. Correlation of m{x) in the case of power-law cor- 
related random mass : Points for each Opu, (power index) are 
almost on one line. 
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FIG. 4. Plots of localization length and typical localization 



length (1): We set L = 300, Op 
the system. Max {|'Ti(a;)|} is about 0.29. 



0.01 and 2000 kinks in 



Next we show the calculations of the typical localiza- 
tion length of the eigenstates. (We obtain "typical" lo- 
calization length by averaging localization lengths of the 
eigenstates within a small range of energy.) As we vary 
the magnitude of the random mass |'7i(x)| in the present 
calculation, we should remark the following facts. If en- 
ergy E is larger than |m(xa)|, the eigenfunction corre- 
sponding that energy is given by the plain wave near 
X ^ Xa- On the other hand if E' < \m{xa)\, the eigen- 
function is a linear combination of exponentially decay- 
ing and increasing functions. Therefore for eigenstate 
with E > Max {|m(a;)|}, the wave function is given by 
a plain wave everywhere and the plain waves are con- 
nected at each position of (anti-)kinks, i.e, jumps of m{x). 
This, of course, does not mean that all eigenstates with 
E > Max {|m(a;)|} are necessarily extended. 

Calculation of the localization lengths vs. energy is 
given for comparatively large in Figs. 4 and 5. There 
are many sharp peaks as in the case of Gaussian dis- 
tributed random mass which was studied in the previ- 
ous paper |^]. To calculate typical localization length 
as a function of energy, we average localization length 
of all states around each peak as in the previous cases. 
The results in Figs. 4 and 5 show the existence of the 
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FIG. 5. Plots of localization length and typical localization 
length (2): We set L = 300, ap„ = 2 and 500 kinks in the 
system. Max {|m(3;)|} is about 0.13. 
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They also show that the above result concerning the 
existence of the mobility edge holds regardless of the 
value of apw ■ Typical localization length is an increasing 
function of energy. This is in sharp contrast to the pre- 
vious result for the system of random kink-distances in 
which the localization length is a decreasing function of 
energy.[] 

We also studied the system size dependence of typical 
localization length. The result is shown in Figs. 6 and 
7. We found that energy of the mobility edge becomes 
larger as we let the system size larger and it can exceed 
the value Max {|TO(a;)|} in large systems. (Max {|m(x)|} 
for different size systems are almost the same.) Local- 
ized states (states whose localization length is smaller 
than the system size) are also observed in the range 
E >Max{|m(a:)|}. 

From the above results, we can conclude that the Max- 
imum value of |TO(a;)| does not give the mobility edge or 
critical energy of the Anderson transition, whereas states 
which have small energy are all localized in large systems. 
(We also expect that the states near E = are always ex- 
tended exceptionally.) This conclusion does not depend 
on the value of am,,- 
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FIG. 6. System size dependence of typical localization 
length (1): We set apw — 1.0. Distance between kinks is 
0.5 and Max {|m(a:;)|} is about 0.10. 
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FIG. 7. System size dependence of typical localization 
length (2): We set Qp„ = 0.1. Distance between kinks is 
0.5 and Max {|m(a;)|} is about 0.20. 



In order to see whether this behaviour of the localiza- 
tion length is due to the long-range correlation of m{x) 
or the way of introducing randomness, we performed the 
following calculations; 

1. Localization length for random values of |m(a;)| 
with exponential correlation and the fixed distances 
between kinks. (We can generate an exponentially- 
correlated m{x) by choosing the modified Bessel 
function as c(m) in Eq.(||) (Fig. 8).) 

2. Localization length for random-kink-distance m{x) 
with fixed |r7i(a;)|. 

The results are shown in Figs. 9 and 10. In the first 
case, the extended states again appear at larger energies 
(see Fig. 9). We see that energy of the mobihty edge 
becomes larger for larger system size. This behaviour is 
the same as in the case of power-law correlated random 
mass. However, the change of the mobility edge seems 
very small. 

On the other hand, (almost) all states are localized in 
the second case (see Fig. 10). This means that the exis- 
tence of the finite mobility edge stems from the quasi- 
periodicity of the random mass with the distances be- 
tween kinks fixed. 

From the above result, we conclude that the appear- 
ance of extended states at high energy is due to the spe- 
cific scheme to generate random mass and it is indepen- 
dent of the correlation of m{x). 



-0.1 I'lf 



0.4 



0.6 



X / System Size 
Correlalion of Random Mass 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
p / System Size 

FIG. 8. An example of m{x) and the correlation: We found 
the exponentially correlated randomness in this example. We 
averaged 500 random systems for the calculation of correlation 
of random mass. 



^It is interesting to study in detail the low-energy behaviour 
of the locahzation length in the present system. 
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FIG. 9. Typical localization length in the system with ex- 
ponentially correlated randomness. Distance between kinks 
is 0.333. Max {|m(a:)|} is about 0.35. 
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FIG. 10. Localization length with exponentially corre- 
lated randomness (with varying kink distances and constant 
|m(a;)|): We set the correlation length A = 1/2, |m(a;)| = 1 
and L — 50. Localization lengths of the states do not exceed 
the system size even in the range E > 1. 



B. Aperiodic random mass 

In the previous subsection, we studied the Dirac 
fermion with random mass of power-low correlation. In 
order to make m{x) correlated algebraicly, we fixed the 
distances between kinks of the telegraphic m(x) and 
made the magnitude of m{x) random. In this subsec- 
tion, we shall consider the second type of random mass 
in which the magnitude of rn(a;) is fixed whereas distances 
between kinks are random variable. 

Recently a closely related model, aperiodic Kronig- 
Penney model in one dimension, was studied in Ref. [|6|. 
The Schrodinger equation of the model is given as 

oo 

- (n/2m)* "(z) + U S{z - Zn)'i'{z) = E •^{z), 

(11) 



n— — oo 



where potential depth U is constant. In the periodic 
case, Zn's (positions of the delta- function potential) are 
arranged periodically, that is, they satisfy Zn — Zn-i = 
A(=constant). We then introduce randomness as the 



displacements of the delta-function potentials and de- 
note them (5„. We define the relative displacement A„ 
as A„ — 5n+i — Sn and also the correlation of A„ as 



^(m) = [AfcA„+„]c„s/[A^] 



(12) 



In Ref. it was shown that extended states (states 
which have vanishing Lyapunov exponent) appear at fi- 
nite energies and therefore the Anderson transition oc- 
curs if £,{m) decreases sufficiently slowly according to the 
power law. We apply this method of generating random 
distances to our model and generate random distances 
between kinks by the following scheme. We first gen- 
erate random variables A„'s from ^(m) by the method 
explained in the previous subsection (we set C =1 and 
k ~ 10^ in Eq.(^) ), and then calculate displacements 
5nS by the following equation, 



(13) 



where C is a parameter which controls the magnitude 
of the randomness. There are three parameters in this 
random system, which are A (mean distance of kinks), 
C (magnitude of the randomness) and decay power of 
^(m) (we denote it v) when we consider the power- law 
correlation. 

In Figs. 11-13 we show C and i(system size) depen- 
dences of the typical localization length. In each case 
the typical localization lengths are a decreasing function 
of E but becomes almost constant as E gets larger (see 
Figs.11-13). 

In the case of large C, the typical localization length 
for nonvanishing E (normalized by the system size L) de- 
creases and become smaller than unity as we let L larger. 
This result indicates that almost all states are localized 
for large C at the infinite system size. We expect that 
this conclusion holds even for small value of u. 



C dependence of Typical Localization Length 
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FIG. 11. C(amplitude of randomness) dependence of typi- 
cal localization length: We set L — 700, A = 0.5 and v = —5.0 
here. In the case of small C, typical localization length at any 
E exceeds the system size. 
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FIG. 12. L(system size) dependence of typical localization 
length: In this figure, we show the L dependence in the case 
of large which means short-range correlated random dis- 
tances. We set A = 0.5, C = 0.05 and u = —5.0 here. 
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IV. CONCLUSION 

In this paper, we studied the Dirac fermion with the 
long-range correlated random mass by using the TMM 
and IVPM. We are especially interested in the delocal- 
ization transition or the existence of a finite mobility 
edge. Recent studies on the Anderson model and Kronig- 
Penney model indicate the existence of finite mobility 
edges. We studied two schemes of the random mass gen- 
eration. In the first one, magnitudes of the mass are 
taken to be random variables with a long-range correla- 
tion whereas distances between kinks are fixed. In this 
case, we found the finite mobility edge. In the second 
one, distances between kinks are random variables. This 
case is closely related with the aperiodic Kronig-Penney 
model in which the finite mobility edge is observed. How- 
ever in the present study, we cannot obtain a definite an- 
swer concerning the (non)existence of the finite mobility 
edge because of the finiteness of the system size. Further 
studies are required for that problem. 

We can also calculate the multi-fractal scaling indices 
directly from wave functions of the random-mass Dirac 
fermion. Values of the indices may depend on the decay 
power of the correlation of the randomness. Results will 
be reported in a future publication |9[ . 



FIG. 13. L dependence of typical localization length: In 
this figure, we show the L dependence in the case of small 
which means long-range correlated random distances. We set 
A = 0.5, C = 0.04 and u = -0.5 here. 

As the value of C becomes smaller, the typical and/or 
mean localization lengths exceed the system size even 
for very large L. Dependence of L is explicitly given in 
Fig. 14. Unfortunately at present, we cannot obtain a 
definite answer to the existence of the mobility edge at 
finite energy. This is in sharp contrast to the Anderson 
model and Kronig-Penney model in which it is believed 
that a genuine Anderson transition occurs. 

L dependence of Averaged Localization Length 
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FIG. 14. L dependence of typical localization length: We 
set A = 0.5 and C — 0.01 here. We show the averaged lo- 
calization length of eigenstates whose energy eigenvalues are 
larger than 0.02. 
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